Determination of energy scales in few-electron double quantum dots 
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The capacitive couplings between gate-defined quantum dots and their gates vary considerably as a function 
of applied gate voltages. The conversion between gate voltages and the relevant energy scales is usually 
performed in a regime of rather symmetric dot-lead tunnel couplings strong enough to allow direct transport 
measurements. Unfortunately this standard procedure fails for weak and possibly asymmetric tunnel cou- 
plings, often the case in realistic devices. We have developed methods to determine the gate voltage to energy 
conversion accurately in the different regimes of dot-lead tunnel couplings and demonstrate strong variations 
of the conversion factors. Our concepts can easily be extended to triple quantum dots or even larger arrays. 



I. INTRODUCTION 

Electrostatically defined coupled quantum dot (QD) 
systems are interesting as an experimental toy model for 
fundamental quantum- mechanical problems f^"^ as qubits 
and registers for quantum informatio n pro cessingPtl^ 
and to simulate molecular electronics. ^^^^ In contrast 
to real molecules, the charge configuration and elec- 
tronic spectrum of these artificial molecules are highly 
tunable p^l^^* the electronic energy scales thus have to be 
redetermined in each experiment. Usually the electronic 
spectra of coupled QDs are measured as a function of gate 
voltages. Hence, a meaningful analysis of the measured 
data involves a conversion of the applied gate voltages 
to energy differences between the electronic states. The 
conversion factors are specific to each sample, and even 
vary if the configuration of a QD system is changed. 

The conversion from gate voltages to energy can be 
achieved by comparison with a known external energy 
scale. The most straightforward method relies on non- 
linear transport measurements where the external energy 
scale is provided by the applied source-drain voltagej^ 
this method w ill be reviewed in Sec. IIVI In the few- 
electron regime ,1^^^ which is desirable for many appli- 
cations, the tunnel barriers of the double QD system 
are often too high to observe a current fiow through the 
double QD. Fortunately, charge fiuct nations can still be 
measured, e.g., using a capacitively coupled quantum 
point contact (QPC) as a charge detector.!^ With this 
technique it is even possible to detect extremely small 
currents through an almost pinched-off double QD in- 
directly by recording tunneling processes in real time.!^ 
This counting method is demanding, though, as it needs 
a special setup, including a low- noise high-bandwidth de- 
tector. A more basic procedure which can be conducted 
with a standard experimental setup would be desirable. 

It is possible to extract energy scales from the thermal 
broadening of the transitions between different charge 
configurations in the stability diagram of coupled QDs^^ 
in an elaborate procedure. This method requires a small 
tunnel coupling to guarantee that the line shapes are 
determined by thermal broadening. In an appropriate 



radio- frequency setup and a double QD with suitable tun- 
nel couplings photon-assisted tunneling can be used for 
energy calibration si nce th e photon energy provides the 
external energy scale.^^^^^ 

Most energy calibration methods including the afore- 
mentioned ones either require special experimental se- 
tups or a narrow regime of tunnel couplings of the QDs. 
In practice the conversion factors are often determined 
once in a rather open double QD by measuring nonlin- 
ear transport, and are then still used after tuning the 
coupled QDs to rather small tunnel couplings. In this 
article we show that such an approach leads to inaccu- 
racies that can be avoided. We present several methods 
to acquire conversion factors, valid for different coupling 
regimes, which are accurate as well as straightforward to 
implement. The discussion here is limited to double QD 
systems with rather weak tunnel coupling between the 
two QDs which, incidentally, have a very large resistance 
compared to the lead resistances; the latter can thus be 
neglected. Strong coupling between the QDs has been 
discussed in Ref. |2l 



II. SAMPLE 

The sample has been fabricated from a GaAs/AlGaAs 
heterostructure containing a two-dimensional electron 
system (2DES) 85 nm below the surface. The charge car- 
rier density of the 2DES was rig = 1.9 x 10^^ cm~^ and 
its mobility was /i = 1.2 x 10^ cm^/(Vs) at low tem- 
peratures. The measurements presented here have been 
performed at an electron temperature of T 50 mK. 
The double QD is defined by applying negative voltages 
to metallic electrodes on the sample surface which have 
been created by electron beam lithography. A scanning 
electron micrograph of a nominally identical structure is 
shown in Fig. ITIa). The design is similar to the one pre- 
sented in Ref.l27l The plunger gates "gL" and "gR" are 
used to change the chemical potentials of the two QDs 
L and R while the tunneling barriers between the QDs 
as well as between each QD and its adjacent lead are 
controlled by four additional gates "tu," "tc," "tl," and 
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FIG. 1. (a) Scanning electron micrograph of a structure 
nominally identical to the actual sample. Eight gates (light 
gray) are used to define a double QD as well as up to two 
QPC charge detectors. Electric circuits for driving current 
through the double QD or the actually used QPC are indi- 
cated, (b) Continuous (black) lines sketch a schematic charge 
stability diagram for an unbiased double QD. The splitting 
of these charging lines, if a source-drain voltage Vdqd is ap- 
plied across the double QD, into source and drain resonances 
is depicted by interrupted (colored) lines (for details see main 
text). The gray triangles mark regions in which first-order 
current is possible. The occupation numbers of the two QDs 
change from (M — 1,N — 1) to M, N from bottom left to top 
right. Overall shifts of the diagram due to the applied bias 
are taken into account assuming equal capacitive couplings 
between the double QD and both leads (see Sec IX for details 
on capacitive coupling), (c) Chemical potentials of the leads 



(/is for source and /xd for drain) and QDs (/xl, Mr) ^he 
positions marked by A, B, and C in (b). (d) Current through 
the double QD measured as a function of plunger gate volt- 
ages V^L and V^R for Vdqd = — 241.1/iV. Regions of finite 
current are enclosed by dashed lines. 



"tr" (these gates also influence the chemical potentials 
but are kept constant during a measurement). Gate "qr" 
together with tr and gR defines the QPC that is used 
as charge detector. A second QPG can be defined on 
the left-hand side of the sample. The voltage V^^r is var- 
ied proportionally to the two plunger gate voltages to 
compensate for the capacitive coupling between plunger 
gates and the QPC and to keep the QPC current roughly 
constant during a measurement. 



III. BASIC RELATIONS 

We define the chemical potential /iQD(A^) of a QD oc- 
cupied by — 1 electrons as the energy needed to add 
the A^th electron. If only a very small bias voltage is ap- 
plied to the double QD, the first-order linear-response 
current, /dqd, through the double QD only flows at 
the triple points where the chemical potentials of both 
QDs (/iL and /ir) and the leads (/is and /xd) are aligned 
(/iL = Mr = Ms = Md)- a charge stability diagram 
containing four different charge configurations (M, A^) is 
shown schematically in Fig. [ijb) in black (continuous) 
lines for negligible bias voltage. Along these charging 
lines the chemical potential of one QD is in resonance 
with the Fermi energy of the leads (/iR = Ms = Md or 

Ml = Ms = Md). 

If a finite bias is applied, resonances with the source 
and the drain lead occur for different plunger gate volt- 
ages. This is also shown exemplarily in Fig. [ijb) where 
we choose Ms > Md- The source resonances Ml,r = Ms 
are plotted as blue (dashed-dotted) lines and the drain 
resonances ml,r = Md as red (dashed) lines. The two 
triple points of finite current (by first-order sequential 
single-electron tunneling) grow into two identical, so- 
called "bias triangles" with /iD S Ml ^ MR ^ Ms, which 
are marked in gray in Fig. [ijb). 1211 The left triangle in Fis. 
[TJb) corresponds to the addition of the first of two elec- 
trons. For its three corners. A, B, and C, the alignment 
of the chemical potentials is sketched in Fig. [TJc). Note 
that the position of the source and drain resonances in 
a stability diagram [such as the schematic in Fig. [ijb)] 
depend directly on /is and md as well as on the capacitive 
couplings between the QDs and the leads (see Sec. [iX] for 
details). 

The chemical potentials mqd of the QDs are con- 
trolled by the voltages Vgi on plunger gates gi via their 
capacitive coupling expressed by the conversion factor 
= — ^QY^ . Considering two plunger gates and a 
double QD, we find 
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From this the slopes 5qd = dVgR/dVgi, of the charging 
lines parallel to AB (CB) in Fig.mb), 



^gR 

^.L 
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as well as the slope of the charge reconfiguration line 
along AC, 



; c^mr = c^ml , 



(3) 
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can be easily derived. In a sufficiently small region of the 
stability diagram all these slopes and conversion factors 
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are constant, which is equivalent to constant capacitive 
couphngs. Under this condition ABC in Fig. [ijb) forms 
a triangle and Eqs. ([l])-([3]) allow an accurate coordinate 
transformation from plunger gate voltages to energy val- 
ues of the QDs' chemical potentials. 



IV. CALIBRATION IN NONLINEAR TRANSPORT 
MEASUREMENTS 

In the simplest case the bias triangles sketched in 
Fig. [ijb) can be measured directly. Such an example 
is depicted in Fig. [TJd) which shows the current flowing 
through the double QD for a bias of Vdqb = — 241.1/iV 
applied to the source lead. Nonzero current is observed 
within two almost identical triangles just as sketched in 
Fig. [ijb) . The conversion factors can be extracted from 
the dimensions of these triangles. We can start, e.g., by 
measuring the voltage change, AI^r, between point B 
and point C along the V^r axis in the stability diagram 
of Fig. [TJd). For the transition B ^ C, we find A/iL = 
and A/iR = eV^QD [compare Fig. [TJc)]. These relations 
and equations ([T])-([3| allow us to determine all relevant 
conversion factors. 
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The numerical values obtained from the current measure- 
ments in Fig. [TJd) are shown in Table |lj Error values are 
also given, which refiect the limited accuracy of determin- 
ing the slopes and voltage differences from the greyscale 
plot. Averaging over both current-carrying triangles has 
been used where applicable, with the error value repre- 
senting the standard deviation. The conversion factors 
presented in Table [T| have to be understood as a lower 
limit of the actual conversion factors due to the limited 
sensitivity of the current measurement, though, as will 



be discussed in detail in Sec. I VIII The error values do 
not include this systematic error. 



V. CALIBRATION IN CHARGE MEASUREMENTS 

Instead of measuring current through the double QD, 
recording its charge by means of a QPC charge detector 
allows a determination of the conversion factors with a 
higher accuracy. In a measurement with zero bias applied 
to the double QD [see Fig.|2|a)] the charging lines (dark) 
and charge reconfiguration lines (white) mark the bound- 
aries between stable charge configurations. For rather 
symmetric tunnel couplings to both leads the triangles 
observed for nonlinear transport may show up again as 
regions of constant average charge^ if a bias voltage is 
applied to the double QD. Such a situation is depicted in 
Fig.|2|b). At finite bias, the white charge reconfiguration 
line is only visible along the short section connecting the 
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FIG. 2. (a) Charge stability diagram (compare Fig. [T]) of 
the double QD measured with the right QPC used as charge 
detector. The transconductance d/gpc/dV^L (determined by 
taking the numerical derivate of the dc current /gpc) is plot- 
ted for (a) Vdqd ^ and (b) Vdqd = -241.1 /xV. Applied 
bias and gate voltages are very similar as in Fig. [ijd). In- 
set: Schematic charge stability diagram similar to Fig. [T]^b). 
Charging lines actually observed in (b) are marked in gray, 
the bias triangles that would be visible in nonlinear transport 
are shaded. 



bases of the two triangles. Additional parallel lines can 
be attributed to excited states in one of the QDs.l^ 



The conversion factors are extracted in the same way 
as for Fig. [l] using Eq. Q and again averaging over both 
triangles. It is helpful to take the slope of the charge 
reconfiguration line from the measurement with zero ap- 
plied bias to increase the accuracy. The resulting slopes 
and conversion factors can be found in Table [H 

If at least one of the three relevant tunnel barriers of 
the double QD (the QD-lead barriers or the interdot bar- 
rier) is increased, /dqd will decrease, eventually becom- 
ing too small to be measured. In this case the energy cal- 
ibration has to be performed via charge detection. The 
features that appear in such a nonlinear stability dia- 
gram depend on the ratio of the three relevant tunnel 
couplings. For instance, the bias triangles are only ob- 
served for rather symmetric tunnel couplings. Fig. [sja) 
depicts a typical situation in which the tunnel coupling to 
the drain lead is small (due to a high barrier) compared to 
the other two relevant tunnel couplings. As a result the 
drain resonance of the left bias triangle is not observed 
and the average charge configuration within this bias tri- 
angle is close to {M^N — 1). The charge reconfiguration 
line (white) is therefore clearly visible in Fig.jsja). From 
the kink in the charging lines, the bias triangles can still 
be reconstructed as demonstrated in Fig. [sj^b), allowing 
the usage of the same calibration relations as described 
above. The results have been added to Table HI 
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SR -0.48 ± 0.03 -0.51 ± 0.03 -0.60 ± 0.04 -0.58 ± 0.03 

SL -1.59 ± 0.18 -1.71 di 0.03 -1.81 =b 0.09 -1.77 ± 0.01 

Scr 1.29 ±0.10 1.48 ±0.12 0.98 ± 0.08 0.70 ± 0.07 



AVgR (mV) 4.52 ±0.08 4.22 ± 0.12 3.86 ± 0.06 5.61 ±0.11 



a^L (meV/V) 


37 ±6 


41 ±6 


56 ±7 


60 ±6 


a% (meV/V) 


76 ±8 


81 ±6 


94 ±7 


103 ±5 


a^j, (meV/V) 


47 ±8 


51 ±6 


53 ±7 


53 ±4 


(meV/V) 


75 ±5 


86 ±7 


96 ±7 


94 ±7 


l^DQD (/iV) 


-241.1 


-241.1 


-241.1 


-387.5 


Vti (V) 


-0.415 


-0.410 


-0.435 


-0.380 


Vtr (V) 


-0.430 


-0.425 


-0.450 


-0.465 



TABLE I. Summary of the charging Une slopes and conversion factors determined from the presented data for the four different 
double QD configurations shown in Figs. [l}j4] Vdqd as well as gate voltages that are changed between measurements are also 
noted for convenience. 
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FIG. 3. (a) Charge stability diagram as in Fig. [2] (Vdqd = 
—241.1 /iV) but with overall reduced and asymmetric tunnel 
couplings. The tunnel barrier between the double QD and 
the drain lead is larger than the other tunnel barriers. Dis- 
tinct kinks in the black charging lines are reminiscent of the 
bias triangles and can be used for determining the conversion 
factors, (b) Same data as in (a); reconstruction of the bias 
triangles (solid lines). 



VI. ASYMMETRICAL CONFIGURATION 



belonging to the transition (M - 1, TV - 1) ^ (M - 1, TV) 
is discontinuous, however, as further away from the triple 
point the cotunneling rate quickly decreases and the right 
QD can be charged more efficiently from the right-hand 
side (source) lead via the large tunnel barrier. The charg- 
ing line of the right QD therefore shows a step as the 
transition from drain to source resonance occurs, with a 
step size determined by e^DgD- In Fig. [4]^c) a sketch of 
the charge stability diagram including source and drain 
resonances is shown, with the lines visible in Fig. 
marked in gray. The chemical potential of the double 
QD system at the points A, B, C, and D are depicted in 
Fig. [4jd). The points C and D have been chosen such 
that the line CD is parallel to the charging line of the 
left QD. At C, the right QD is in resonance with the left 
lead whereas at D the right lead is resonant, so A/iL = 
and A/iR = eV^QB between C and D, just as between 
B and C in Figs. [l]-[3] Hence, the distance CD deter- 
mines AVgYi' To obtain the conversion factors listed in 
Table [l| from Eq. Q we additionally use the slopes of the 
charging lines, 5l, sr, and Scr- 



Finally we study an even more asymmetric double QD 
system, with the barrier between the right QD and the 
adjacent source lead almost closed. Such a situation is 
favorable for experiments studying, e.g., backaction of 
the QPC on the double QD^^^. Current through the 
double QD is almost completely blocked (/dqd — 0). 

Fig. |4ja) shows the charge stability diagram of such a 
double QD system with zero bias voltage in comparison 
to Fig. ^h) where a large negative bias is applied to the 
right-hand side (source) lead. No indications of triangles 
can be observed in Fig.EJb), though. Most of the charg- 
ing lines visible in Fig. HTa) are resonances with the left 
(drain) lead since the left tunnel barrier is reasonably low. 
In the vicinity of the triple points (//l — Mr) charging of 
the right QD occurs from the left (drain) lead in an elas- 
tic cotunneling process via the left QD. The charging line 



VII. COMPARISON OF RESULTS 

The conversion factors listed in Table [T| vary consid- 
erably between the four measurements, even though the 
data shown in Figs. [l]-|4] have all been measured within 
the same cooling run on the same sample. The main 
difference between those measurements is that the dou- 
ble QD is tuned to slightly different geometries by vary- 
ing the gate voltages Vn and Vtr as shown in Table [T| 
(Vtu = -0.440 V and Vtc = -0.310 V are kept constant). 
For instance the conversion factors of Fig. [4] are up to 
62 % higher than those of Fig. [T] We conclude that the 
capacitances between gates and the QDs depend on the 
detailed QD geometries which are strongly affected by the 
applied gate voltages. Whenever accurate energy values 
of the electronic states in coupled QDs are desired, it is 
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FIG. 4. (a) Charge stability diagram of a double QD with 
very asymmetric tunnel couplings at Vdqd — 0. The tunnel 
barrier between double QD and source lead is almost closed, 
(b) Measurement as in (a) but for Vdqd = —387.5 /iV. (c) 
Schematic charge stability diagram as in Fig. mb). Charging 
lines are underlined in gray where visible in (b). The volt- 
age difference AVgR between positions C and D can be used 
for determining the conversion factors (for details see main 
text), (d) Chemical potentials for the positions marked A, 
B, C, D in (c). The source resonance in "A" is not observed 
in (b), since charge exchange through the very large source 
barrier (denoted by the thick black line in the level diagrams) 
is almost blocked. 



therefore of utmost importance to perform a calibration 
right at the gate voltage settings of interest. In this arti- 
cle we demonstrate that simple calibration methods are 
available for very different coupling regimes. 

The voltage settings for Figs, [l] and [2] are almost iden- 
tical, and the same should be true for the conversion fac- 
tors. Nevertheless the a's listed in Table [l| are up to 14 % 
higher for Fig. [2] than those determined from Fig.[l] This 
discrepancy is caused by the difficulty and according in- 
accuracy in determining the size of the current-carrying 
triangle in Fig. [l] since the observed size depends on the 
sensitivity of the current measurement. This causes a 
systematic error in AVgR which is not included in the 
error value given in Table [T[ 

In addition, in a current measurement the size of the 
bias triangle is influenced by the width of the transi- 
tions between the different charge configurations (i.e., the 
width of the charging lines in a charge measurements) 
which depends on the tunnel couplings. The accuracy 



is much higher for charge detection, where AV^r can be 
determined by using the center of the charging lines in- 
dependent of their width. Hence, we recommend using 
charge detection or, at least, the derivative dl^QB/dYgi 
of the current for energy calibration. 



VIII. INTERPLAY OF THE DIFFERENT METHODS 

In Fig. [5] a larger portion of a charge stability dia- 
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FIG. 5. Larger section of the charge stability diagram for 
an asymmetrically tunnel coupled double QD similar to that 
of Fig. [I^b) at yDQD = -387.5 /xV. Absolute values of the 
electron numbers in the left (L) and right QD (R) are given 
as L/R. Barely visible almost horizontal lines are related to 
resonances in the detector QPC. 

gram (measured using charge detection) is shown which 
features all of the effects described in Sections IVT^VII as 
the ratio between the relevant tunnel couplings varies 
throughout the stability diagram as a function of the 
plunger gate voltages. For this measurement the dou- 
ble QD has been tuned rather asymmetrically, similar 
to the configuration in Fig. |4j and the feature visible in 
the latter graph reoccurs here as the electron numbers 
change from 1/1 to 2/2. As the tunnel couplings vary, 
different effects are observed at other transitions. The 
feature visible at the transition 0/0 ^ 1/1 is very similar 
to that of Fig. [2] while the transition 0/1 ^ 1/2 resem- 
bles Fig. [sj The features visible at 1/0 ^ 2/1 contains 
mixed signatures. The region 2/2 ^ 3/3 shows no signs 
of ^DQD 7^ 0. Here our calibration methods fail and the 
conversion factors have to be estimated by extrapolation 
from other areas of the stability diagram. Calibration 
methods which need specific experimental setups might 
still work (compare Sec.[l|). 
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IX. CAPACITIVE COUPLING BETWEEN QUANTUM 
DOTS AND LEADS 

All calibration methods described so far rely on dis- 
tinct features that appear in the stability diagram as a 
finite but fixed bias Vbqb is applied. It is tempting to 
utilize instead the linear response of the position of a 
charging line as Vbqb is varied. Fig. [6] illustrates the 
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FIG. 6. The gray-scale plot shows the zero-bias data from 
Fig. |4ja), with the charging lines from Fig. (Vdqd = 

—387.5 /iV) superimposed as dashed lines. 

overall effect of applying a bias voltage. The zero-bias 
data from Fig. |4]^a) are reproduced, with the charging 
lines from the Vbqb = —387.5 /i data in Fig. |4|b) super- 
imposed (dashed lines). Changing the bias results in a 
much smaller shift of the charging lines than expected 
from the split charging line on the left-hand side that 
has been used for calibration in Sec. |Vll The observed 
shift is therefore not a direct measure of eV^QB- It is 
much smaller than expected from the applied bias Vbqb 
because of a compensation related to the electrostatic 
coupling between the QDs, the gates, and the leads. As 
pointed out with this example, relying solely on changes 
in the charging line positions is not suitable to obtain 
a valid calibration. Instead, we used features that show 
transitions from source to drain resonances. 



X. CONCLUSIONS 

We have described several methods of obtaining the 
conversion factors between gate voltages and chemical 
potentials in a double QD system in the few-electron 
regime. Depending on the specific couplings and tunnel- 
ing rates in the system, different phenomena are observed 
in the charge stability diagram as a bias is applied across 
the double QD, and most of them can be harnessed to 
determine the conversion factors. Even for a very asym- 
metric system, a calibration procedure has been devel- 
oped. The latter method might gain significance if the 
number of quantum dot in series is increased, as addi- 
tional QDs tend to act in the same way as large barriers 
when in Coulomb blockade. The procedure might there- 
fore be extended to few-electron triple QD circuits which 
include a charge detectoi^^ 
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